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Abstract 

Spin-density-functional theory of quantum point contacts (QPCs) reveals the 
formation of a local moment with a net of one electron spin in the vicinity of 
the point contact - supporting the recent report of a Kondo effect in a QPC. 
The hybridization of the local moment to the leads decreases as the QPC 
becomes longer, while the onsite Coulomb-interaction energy remains almost 
constant. 
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The discovery of the extra conductance plateaus at 0.7(2e 2 //i) in quantum point con- 
tacts (QPCs) and at 0.5(2e 2 //i) in clean quantum wires [|J, in addition to the usual 
quantization of the conductance into steps of 2e 2 /h has focused attention on the role 
of electron-electron interaction in these low- dimensional quantum systems |7|,p|,P, p!0|JTl| , p!2[] . 

A recent experiment on QPCs has revealed several features characteristic of the Kondo 
effect [|13| . Specifically, there is a zero-bias peak in the differential conductance which splits 
in an in-plane magnetic field |14|]. A single energy scale, the "Kondo" temperature Tk, 
sets the width of the zero-bias peak, the magnetic field required to split the peak, and the 
crossover temperature to perfect transmission. Moreover, the extra conductance plateau 
at 0.7(2e 2 //i) can be explained using an Anderson model in which the hybridization of a 



localized electron to the band depends on energy and valence fl5|| . A puzzling question is 
how a localized spin can form in an open QPC system. 

In this work, we employ spin-density-functional theory (SDFT), and find that a local 
moment with a net of one electron spin is formed in the vicinity of a QPC. The splitting 
between the two spin directions at high in-plane magnetic fields JTJJ approaches a finite 
residual splitting at zero field, which we interpret as the Coulomb-interaction energy U 
between electrons at the site of the QPC. ^From the width of the localized state we infer 
that the hybridization to the leads decreases for longer QPCs, while the Coulomb U remains 
nearly constant. Our interpretation that a dynamical local moment forms at the QPC differs 
from that of Wang and Berggren [0 who infer a true spin-polarized ground state. We believe 
that the experimental observation of signatures of the Kondo effect Jl3|], combined with a 



theorem forbidding spin-polarization in one dimension fll6| , strongly favor the dynamical- 
spin interpretation. 

The electronic states of a realistic quantum point contact are obtained using density- 
functional theory within the local-density approximation ||17| . This method allows us to 
study both potential confinement and electron-electron interaction in a unified framework. 
First, for a clean quantum wire uniform in the x-direction and with a parabolic confining 
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potential in the y-direction of V^ d (y) = {l/2)m*uj 2 y 2 \ we readily find the wave functions 
e ±tk ™ x i/}°(y), chemical potential /i, band-edge energies e^, and the charge density p 0a (y) for 
a given density n\D [|18j. Here n is the subband index and a is spin index, respectively. 
Next, we introduce the QPC potential, 



V QPC (x, V) = V(x)/2 + m* (uj y + V{x)/hf y 2 /2 - V&fo) 



(1) 



(see insets to Fig. 1), where V(x) = Vo/cosh 2 (x/<i), with decay length d = J2Vo/m* /u x , 
and solve the following Kohn-Sham equation for wave functions with energy e, 



h 2 d 2 



2m* dx 2 
+ gfiBB\\<T + e° 



+ V QPC {x, y) + SV H (x, y) + 6V£(x, y) 

Kk n ( x >y) = e Kkn( x >y)- 



(2) 



Electrons incident from the two leads (x — > ±oo) are scattered elastically by the effec- 
tive QPC potential 5V a (x, y) = Vqpc(x,y) + 8Vh(x, y) + 8V£ c (x, y), which is the difference 
between the QPC self-consistent potential and that of the clean wire. The Hartree and 
exchange-correlation parts of the potentials are written, respectively, as 

„2 



5V H (x,y) 



dx'dy'{p(x',y')-p°(y f )} 
1 



(x — x') 2 + (y — y') 2 y( x ~ x ') 2 + (y ~ V 1 ) 2 + a<2 
6E xc [p,^] 6E xc [p°,e] 



6V: c (x,y) 



(3) 
(4) 



5p a (r) 5p 0a (y) ' 

where we have introduced an image-charge plane at a distance a = lOOnm to model the 
experimental geometry | I8| . The exchange-correlation energy functional E xc [p,^] is treated 
in the local density approximation with the local spin polarization £(r) = [p (r) — p^(r)]/p(r) 



The eigenstates of (fj) can be characterized as waves incident from the left k a(x,y) 
and from the right $?* k Jx,y). Expanding, ^ k Jx,y) = E m ^, ra ,te(s)C(!/)i 



the form of plane waves with wavevectors k a n = J2m*(e — €°)/H far from the QPC region: 
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U n,m,kZ \ X ) 



fD ik n x £ _j_ ^tj p ik m x . ™ <^ — ™ 
c u n,m i 1 n rnF > — 

( 5 ) 



where m and m are ^ ne elements of unknown reflection and transmission matrices, and 
where \xo\ = 500nm is sufficiently far from the QPC that SV(x,y) is negligible. ^/^ k <y(x,y) 
is expanded analogously. We employ the recursion-transfer-matrix (RTM) method [21J to 



solve the Kohn-Sham equation (|2[) with the above boundary conditions. ^From the resulting 
wave functions, the density p(x, y) is constructed as a sum over occupied states, 



1 r°° ( f> h a2 \ tit, 

L/R 



(6) 



where the Fermi distribution function is /(e) = l/[exp {(e — fi)/kBT} + l]. These procedures 
are iterated until self-consistent solutions are obtained for the electron density nuj. We use 
the material constants for GaAs, m* = 0.067mo, n = 12.9, g = 0.44 and fix the external 
confinement in the y-direction in the wire such as %u y = 2.0meV. 

Figure 1 shows the electronic properties of QPCs with Vq = 3.0meV at T = 0.1K for three 
different lengths: (a),(d) d = 82.6nm (fru; x .=1.0meV), (b),(e) d = 55.0nm (^.w^^l.SmeV), 
and (c),(f) d = 41.3nm (hu x =2.0meV). The density far from the QPC is taken as n\D = 
2.80 x 10 _2 nm _1 , where the electrons far into the wire are unpolarized and only the lowest 



two spin subbands contribute to transport |18 |. We find that a solution with broken spin 



symmetry coexists with an unpolarized solution. To obtain the broken-symmetry solution, 



we first apply an in-plane magnetic field of up to B = 6T |T4]] , and then solve the Kohn-Sham 
equations self-consistently while reducing the magnetic field to zero. 

In Fig. l(a)-(c) we show the self-consistent QPC barrier as a function of position x in 
the direction of current flow. Specifically, we plot the energy of the bottom of the lowest 
ID subband Cq(x) relative to the band-edge far into the wire, for both spin- up (solid 
lines) and spin-down (dashed lines) electrons. Small Friedel oscillations with a period of 
27r/2/c/ ~ 72nm are present far into the wire. The self-consistent QPC barrier is strongly 
spin dependent in all three cases, with the chemical potential /i lying above the spin-up 
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barrier but below the spin-down barrier. These observations are consistent with the results 
of Wang and Berggren p| . 

In Fig. l(d)-(f) we show the ID electron density in the QPC. The solid lines give the 
net spin-up density and the dashed lines give the spin-averaged densities, which approach 
uid/2 = 1.40 x 10~ 2 nm _1 far into the wire. For all three QPC lengths, there is an excess 
spin-up density in the vicinity of the barrier, with Friedel oscillations persisting into the 
wire. In each case, the integrated net spin-up density is close to one spin: 0.85 for (a), 0.93 
for (b), and 0.90 for (c). Thus there is a local moment with a spin of 1/2 formed at the 
QPC. 

The local moment found within SDFT results from the self-consistent flattening of the 
QPC barrier. For example, above a square barrier there is a series of quasi-bound states, 
resulting from multiple reflections from the edges of the barrier. While the bare QPC 
potential does not give rise to resonances, the self-consistent potential is flattened for spin- 
up, particularly for the longer QPCs, and the first resonance in the spin-up local density of 
states is clearly resolved. Since this resonance lies below the chemical potential, it is fully 
occupied. The result is a localized net spin of 1/2 at the QPC. 

An important question is whether SDFT is reliable in predicting a local moment. SDFT 
is a mean-field theory, and has limited utility for strongly correlated electron systems. In 
particular, consider an Anderson model consisting of a single site with a Coulomb interaction 
U between two spin-degenerate orbitals of energy Eq, with hybridization T to band electrons 
|22|]. For partially occupied orbitals, the local density of states for each spin splits into two 



peaks one at Eq and the other at Eq + U (and at low temperatures an additional Kondo peak 
at the chemical potential) [p3[| . Taken at face value, SDFT is inadequate for this model as 
it can give only one peak for each spin. However, SDFT has been applied successfully to 
calculate local moments in a system closely corresponding to the Anderson model, namely 
transition-metal adatoms on surfaces ||24| . The key is that even though the SDFT solution 
breaks spin-rotation symmetry, the frozen magnetization can still be reliably interpreted 
as the magnitude of the dynamic local moment. In practice, we follow the SDFT method 



developed by Janak to find the magnetization of metallic metals [ 25| . Specifically, we first 



solve the SDFT equations in a polarizing magnetic field, and then reduce the field to zero 
to extract the parameters U and V of the underlying Anderson model. 

To estimate the properties of the bound state, we go beyond the formal validity of SDFT 
and study the Kohn-Sham orbitals. Specifically, we plot the local ID density of states v a {e) 
at the center of the QPC in the insets of Fig. l(a-c), for the respective physical parameters. 
Below the chemical potential, the local density of states ^t(e) for spin- up electrons shows a 
resonance which broadens as the QPC is shortened. 

Figure 2 depicts the local ID density of states ^o-(e) for huj x = 1.5meV in the presence 



of in-plane magnetic fields B from OT to 10T in steps of IT pHfl . Traces are vertically 
offset by a constant amount. The solid lines are for spin- up and the dashed lines are for 
spin-down electrons. With increasing Zeeman splitting, the resonance for spin-up electrons 
shifts to lower energies while the onset of the spin-down density of states shifts to higher 
energies. From the residual splitting between these features at B = we obtain an estimate 
for the Coulomb-interaction energy U in a local-moment description of the QPC. Similarly, 
we obtain the hybridization V from the width of the spin-up resonance. 

Fig. 3 shows the dependence of the Coulomb energy U and the hybridization T on the 
length of the QPC. U is obtained from the energy difference between the resonance center of 
z^(e) and the energy at which the derivative of v±(e) is a maximum. T is obtained from the 
FWHM of a Lorentzian fit to the resonance in ^(e). The hybridization energy T increases 
sharply from T ~ O.lmeV up to ~ 0.6meV as the QPC is shortened from d = 80nm to 
d = 40nm. The Coulomb energy U, on the other hand, stays nearly constant at ~ 0.6meV 
over this range. 

Within the Anderson model, a strong Kondo peak in the differential conductance is ex- 
pected when kT < ksT K < Y <C U. This suggests that there is an optimal range of QPC 
lengths for observing the Kondo effect: for very long QPCs, T and therefore the Kondo 
temperature Tk will be too small compared to the real temperature, whereas for very short 
QPCs, the resonance width T will be too broad for a Kondo effect ever to develop. Exper- 
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imentally, Reilly et al. observed that the extra conductance plateau in QPCs decreased 
from 0.7(2e 2 //i) to 0.5(2e 2 //i) with increasing QPC length, suggesting a suppression of the 
Kondo effect [T5[]. Our SDFT results suggest that it is primarily the decreasing hybridiza- 
tion of the local moment to the leads that is responsible for the decrease of the Kondo 
temperature and the evolution of the conductance plateau. 

We have presented SDFT results only for chemical potentials sufficiently above the res- 
onance energy that the local moment is fully formed. SDFT is generally unreliable for a 
partially filled orbital, unless self-interaction effects are explicitly removed [P6|| , and this is 
not practical for an open system with long range interactions like a QPC. 

In conclusion, we studied the electronic states of quantum point contacts (QPCs) using 
the spin-density-functional method. We found that a local moment with a spin of 1/2 is 
formed in the vicinity of the QPC barrier. This strongly supports recent claims of a Kondo 
effect in transport through a quantum point contact ||13|| . For the local moment, we obtained 
estimates for both the onsite Coulomb-interaction energy and the hybridization to the leads. 
The latter decreases rapidly with increasing length of the QPC. 
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DOE grant DE-FG02-00ER45790. 
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Figure Captions 

Fig.l: Results for QPC potentials of increasing sharpness: (a),(d) huo x = l.OmeV, (b),(e) 
huj x = 1.5meV, and (c),(f) %ui x = 2.0meV. (a)-(c) Self-consistent "barrier", i.e. energy 
of the bottom of the lowest ID subband e (a;) at temperature T = 0.1K as a function 
of position x in the direction of current flow through QPC. The chemical potential 
li is indicated by an arrow on the left. Solid lines are for spin- up and dashed lines 
are for spin-down electrons. Insets: local density of states v a {e) at the center of the 
QPC. (d)-(f) ID electron density in QPC. The solid line gives the net spin-up density 
and the dashed line gives the spin-averaged density. Insets: contour plot of the QPC 
potential V Q p C (x,y). 

Fig. 2: Local density of states at the center of the QPC for in-plane magnetic fields B\\ from 
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T to 10 T in steps of 1 T [14]. Traces are vertically offset by a constant amount. 



The solid lines are for spin-up and the dashed lines are for spin-down electrons. 

.3: Hybridization energy T and on-site Coulomb energy (7asa function of hu x , the sharp- 
ness of the external QPC potential. V is obtained from the FWHM of the Lorentzian fit 
to the resonance. U is obtained from the energy difference between the resonance cen- 
ter and the maximum of du(e)/de for the high-energy spin. Inset: ^(e) and dv^{e)/de 
at the center of QPC for Tioj x = 1.5meV. The dashed line is the Lorenzian fit to 
the resonance of ^-f(e). Hybridization energy T and on-site Coulomb energy U are 
indicated. 
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